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This  study  was  undertaken  to  investigate  the  implication  of  geoacoustic  behaviors  in 
the  shallow  marine  sediments  associated  with  the  changes  in  geotechnical  index 
properties.  Two  piston  cores  (270  cm  and  400  cm  in  core  length)  used  in  this  study 
were  recovered  from  stations  1  and  2,  the  western  continental  margin,  the  East  Sea. 
Scanning  electron  microscopy  (SEM)  was  employed  to  illustrate  the  effects  of 
microstructure  on  shear  properties.  The  direct  SEM  observation  of  sediment  fabrics 
is  inevitable  to  understand  the  correlation  of  the  changes  in  geoacoustic  properties  to 
the  sediment  structure.  The  consolidation  of  sediments  by  overburden  stress  resulting 
in  the  clay  fabric  alteration  appears  to  play  an  important  role  in  changing 
shear  properties.  Water  contents  and  porosity  of  sediments  gradually  decreases  with 
increasing  depth,  whereas  wet  bulk  density  shows  a  reverse  trend.  It  is  interesting  to 
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note  that  shear  wave  velocities  increase  rapidly  from  8  to  20  m/ s  while 
compressional  wave  velocities  significantly  fluctuate,  ranging  from  1450  to 
1550m/s  with  depth.  The  fabric  changes  in  sediment  with  increasing  depth  for 
example,  uniform  grain  size  and  well  oriented  clay  fabrics  may  cause  the  shear 
strength  increase  from  1  to  12kPa.  Shear  wave  velocity  is,  therefore,  shown  to  be 
very  sensitive  to  the  changes  in  undrained  strength  for  unconsolidated  marine  sedi¬ 
ments.  This  correlation  allows  an  in-situ  estimation  of  shear  stress  in  the  subsurface 
from  shear  wave  velocity  data. 

Keywords  geotechnical  and  shear  properties,  microstructure,  SEM,  the  East  Sea 


For  the  past  three  decades  geotechnical  and  geoacoustic  properties  of  marine  sedi¬ 
ments  have  been  of  a  great  interest  to  various  fields  such  as  geophysics,  seafloor 
engineering,  sedimentology,  soil  mechanics,  and  underwater  acoustics  (Hamilton 
1971,  1980;  Hamilton  and  Bachman  1982;  Richardson  et  al.  1991,  1997,  2002). 
The  history  of  sea  level  changes  can  be  revealed  by  the  information  of  textural 
changes  and  the  burial  diagenetic  deformation  obtained  by  acoustic  properties 
(Hamilton  1980;  Hamilton  and  Bachman  1982).  Shear  wave  velocity,  shear  modu¬ 
lus,  and  shear  strength  are  used  to  predict  the  stability  of  sediment  slopes,  consoli¬ 
dation  behavior  of  sediments,  and  the  strength  of  marine  foundation  (Nacci  et  al. 
1974;  Theilen  and  Pecher  1991).  This  is  because  the  shear  wave  velocity  in  uncon¬ 
solidated  marine  sediments  has  been  shown  to  be  sensitive  to  overburden  pressure 
or  compaction  (Brunson  1993).  The  stability  of  seafloor  sediments  is  closely 
related  to  the  physical  state,  behavior  to  the  external  loading,  and  the  stress  his¬ 
tory  (Davis  et  al.  1993).  The  depositional  environment  and  sediment  supply 
initially  control  the  texture,  sediment  grade  and  sedimentation  rate.  These  charac¬ 
teristics  in  turn  influence  the  petrophysical  properties  such  as  porosity,  water  con¬ 
tent,  and  grain  packing.  The  consolidation  of  sediments  has  been  causing  a 
porosity  reduction  by  mechanical  overburden  pressure  with  increasing  depth, 
which  is  well  documented  (Kim  et  al.  1999,  2001;  Bryant  et  al.  1991).  The  weight 
of  the  overlying  sediments  rearranges  the  individual  sediment  particles  (fine  clays) 
and  domains  manipulating  the  petrophysical  properties  of  sediments  (Bennett  et  al. 
1991;  Bryant  et  al.  1991). 

Previous  studies  have  given  little  consideration  on  the  relation  between  sediment 
fabric  and  shear  properties,  which  can  be  used  to  develop  a  model  to  predict  an  in- 
situ  estimation  of  shear  stress.  Most  studies  have  mainly  focused  on  the  change  of 
microstructure  related  to  geotechnical  index  properties  (Bowles  et  al.  1969;  Bohlke 
and  Bennett  1980;  Bennett  et  al.  1981,  1989;  Kim  et  al.  2001),  and  little  is  known 
about  the  effect  of  microstructure  on  shear  properties,  particularly  of  shear  wave 
velocity  in  marine  sediments.  In  this  article,  the  response  of  shear  wave  velocity  to 
fabric  changes  is  studied  by  direct  observation  of  the  particle  arrangements  of  sea¬ 
floor  sediments.  This  suggests  that  the  relationship  between  shear  properties  and 
microstructure  is  important  to  understand  the  geoengineering  behavior  of  marine 
sediments. 

Materials  and  Methods 

Two  piston  core  samples  used  in  this  study  were  recovered  from  station  1  (Water 
depth:  2,205  m,  Lat.:  36°52’  6248N,  Long.:  130°07,547E)  and  station  2  (Water  depth: 
1,481  m,  Lat.:  36°06TN,  Long.:  130°05,24E),  the  western  continental  margin  of  the 
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Figure  2.  Columnar  sections  of  the  cores  at  station  1  and  2.  These  cores  consist  of  homo¬ 
geneous  mud  with  tephra  layer  (St.  1)  and  bioturbation  (St.  2).  Mean  grain  size  is  also  dis¬ 
played.  Arrow  marks  of  station  1  are  the  depth  selected  for  SEM  micrograph. 


The  wet  mass  (g),  dry  mass  (g),  wet  volume  (cm3)  and  dry  volume  (cm3)  of  the 
sample  were  measured  using  methods  described  by  Blum  (1997).  The  samples  were 
oven-dried  at  105°C  for  24  hours  and  allowed  to  cool  in  a  desiccator  for  about 
3  hours.  Wet  and  dry  volume  of  samples  were  determined  using  a  Quantachrome 
Penta-Pycnometer,  helium-displacement  pycnometer.  Sample  volume  was  determ¬ 
ined  at  least  five  times,  until  readings  were  consistent  (i.e.,  standard  deviation  less 
than  0.01%).  Sample  mass  was  determined  with  an  error  within  0.1%  using  Scitech 
electronic  balances. 
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Geoacoustic  Properties 

Compressional  and  shear  wave  velocity  was  measured  on  the  split  cores  from  top  to 
the  bottom  at  10-cm  interval  using  a  pulse  transmission  technique  (Birch  1960).  To 
improve  the  coupling  between  the  transducer  head  and  the  sample,  distilled  water 
was  applied  to  the  contact.  Driving  frequency  and  voltage  were  set  at  75  kHz  and 
5  V  p-p  for  P-wave;  500  Hz  and  10  V  p-p  for  S-wave,  respectively.  Pulse  generator 
(Model:  Wavetek  178,  50  MHz)  and  oscilloscope  (Model:  LeCroy  9400A,  Dual 
175  MHz)  were  used  to  generate  signals  and  detect  the  signals  transmitted  through 
the  sample.  A  bimorph  ceramic  bender  transmitter  and  receiver,  designed  and  con¬ 
structed  at  the  Naval  Research  Laboratory-Stennis  Space  Center,  were  utilized  to 
measure  shear  wave  velocity.  Shear  wave  velocity  was  measured  along  both  the  hori¬ 
zontal  and  the  vertical  directions  to  the  core-axis. 


Scanning  Electron  Microscopy  (SEM) 

The  initial  sample  pretreatments  were  performed  following  the  L.R.  White  (LRW) 
resin  impregnation  method  described  by  Kim  and  Peacor  (1995).  This  method  allows 
the  original  fabric  of  the  sample  to  remain  intact.  The  samples  (50  cm  and  214  cm  in 
core  depth  from  station  1)  were  polished  by  using  diamond  pastes  and  then  coated 
with  gold.  The  significant  variations  in  shear  strength  and  shear  velocity  were 
observed  in  these  depth  ranges  of  station  1 .  SEM  observations  were  carried  out  with 
a  HITACHI  S-2400  at  20  kV  and  a  beam  current  of  80-90  nA  in  the  cooperative  lab¬ 
oratory  center  at  Pukyong  National  University,  Busan,  Korea. 


Results 

Grain  Size  Analysis 

The  sediment  in  station  1  shows  a  coarse-grained  layer  (about  6.50),  tephra  layer  ori¬ 
ginated  from  the  Ulleung-Oki  volcanic  eruption  (9.3  ka  ago)  at  the  210  cm  in  depth 
(Figures  2  and  3).  Gravel-sized  pumices  also  occurred  at  this  interval,  while  silt-sized 
particles  are  shown  around  the  62-cm  depth  (Figures  2  and  3).  The  sand  and  silt¬ 
sized  particles  (e.g.,  pumice,  foraminifers)  are  scattered  around  the  230-cm  depth 
and  the  mean  grain  size  slightly  increases  from  8.5  to  80  with  increasing  depth. 
The  average  contents  of  sand,  silt,  and  clay  are  1.79%,  24.51%,  72.27%,  respect¬ 
ively.  The  sediment  in  Station  2  is  characterized  by  the  uniformed  grain  size  (about 
8.80)  and  homogeneous  mud  with  bioturbations  (Figures  2  and  4).  The  average  con¬ 
tents  of  sand,  silt,  and  clay  are  0.06%,  20.47%,  79.47%,  respectively. 


Geotechnical  Index  Properties 

Water  content  and  porosity  of  the  sediment  in  station  1  show  an  abrupt  decrease  at 
62  cm  (from  65  to  50%  in  water  content,  from  85  to  75%  in  porosity)  and  230  cm 
(from  60  to  50%  in  water  content,  from  80  to  70%  in  porosity)  in  depth  (Figure  3), 
which  may  be  caused  by  particle  size  increase  (coarse  grains),  while  wet  density  shows 
a  reverse  trend  of  water  content  and  porosity.  Dry  density  shows  the  lowest  value 
(about  2.1g/cm3)  at  210cm  depth  (Figure  3)  where  a  tephra  layer  (e.g.,  pumice) 


42 


G.  Y.  Kim  et  al. 


Mean  grain  size  (o) 
6  7  8  9  10 


Figure  3.  Mean  grain  size,  geotechnical  properties,  and  p-wave  velocity  for  core  samples  at 
station  1.  The  values  are  markedly  changed  the  intervals  below  210  cm  depth. 

appears  and  the  highest  value  (2.57  g/cm3)  around  the  240-cm  depth  due  to  the  abun¬ 
dant  sand-sized  foraminifers. 

Water  content  and  porosity  of  the  sediments  in  station  2  show  a  similar  trend 
with  each  showing  approximately  60%  and  80%,  respectively  (Figure  4).  As  shown 
in  Figure  4,  the  highest  values  appear  at  80  cm  depth.  Wet  density  shows  a  mirror 
image  of  water  content  variation  having  the  lowest  value  at  80  cm  depth  correspond¬ 
ing  to  a  higher  water  content  and  porosity.  The  variation  of  dry  density  shows  from 
2.2  to  2.5  g/cm3  with  depth. 

Geoacoustic  and  Shear  Properties 

Compressional  wave  velocity  measured  in  sediments  of  station  1  shows  fluctuations 
from  1450  to  1550  m/s  (Figure  3)  and  shear  wave  velocity  is  significantly  increased 
from  approximately  5  m/s  to  14  m/s  at  the  depth  of  100  cm  (Figure  5).  There  is  no 
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Figure  4.  Mean  grain  size,  geotechnical  properties,  and  p-wave  velocity  for  core  samples  at 
station  2.  The  values  are  characterized  by  an  uniform  variation  with  depth. 
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Shear  wave  velocity  (m/s) 
5  10  15  20  25  30 


Void  ratio  Shear  strength  (kPa) 
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Figure  5.  Shear  wave  velocity,  void  ratio,  and  shear  strength  versus  depth  at  station  1 .  The 
values  are  significantly  increased  from  the  boundary  of  100cm  depth. 

systematic  variation  of  compressional  wave  velocity  for  the  core  from  Station  2,  and 
the  pattern  is  also  similar  to  that  in  Figure  3. 

Shear  properties  (shear  wave  velocities,  shear  strength)  and  void  ratio  gradient 
are  plotted  with  depth  (Figures  5  and  6).  Shear  wave  velocities  in  the  horizontal 
and  vertical  directions  show  a  similar  pattern  with  depth.  Shear  wave  velocity  of 


Shear  wave  velocity  (m/s) 
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Figure  6.  Shear  wave  velocity,  void  ratio,  and  shear  strength  versus  depth  at  station  2.  The 
values  gradually  increase  with  depth. 
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Figure  7.  Clay  fabric  of  sample  taken  at  a  depth  of  50  cm  at  station  1.  Top  (a):  SEM  micro¬ 
graph  showing  typical  clay  domains  and  clay  particles.  Middle  (b):  SEM  micrograph  showing 
framboidal  pyrite  cluster,  skeletal  grains,  and  clay  particles.  Bottom  (c):  EDS  spectrum  show¬ 
ing  chemical  composition  (FeS2)  of  framboidal  pyrite.  Note  that  iron  and  sulfide  are  greatly 
enriched. 
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Figure  8.  Clay  fabric  of  sample  taken  at  a  depth  of  21 5  cm  at  station  1 .  Top  (a)  and  middle  (b):  SEM 
micrographs  showing  parallel  orientation  and  well  oriented  clay  fabric.  Bottom  (c):  SEM  micro¬ 
graph  showing  typical  cardhouse  clay  fabric,  the  voids  appear  at  the  center  and  upper  right. 
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(approximately  10  m/s)  than  that  of  the  model  (i.e.,  shear  wave  velocity  gradient 
with  sediment  depth)  suggested  by  Richardson  et  al.  (1991).  Hamilton’s  (1987) 
and  Bryan  and  Stoll’s  (1988)  models  have  higher  values  than  those  of  Richardson 
et  al.  (1991)  for  a  similar  sediment  type.  This  discrepancy  may  be  caused  by  environ¬ 
mental  parameters,  such  as  sedimentary  processes  and  environments,  various  sedi¬ 
ment  types,  degree  of  consolidation  and  compaction,  and  burial  depth 
(Richardson  et  al.,  1991). 

The  changes  in  shear  wave  velocity  caused  by  increased  overburden  pressure  are 
noticed  in  Figures  5  and  6.  The  value  was  significantly  increased  from  8  to  approxi¬ 
mately  25  m/s  (Figure  6),  suggesting  that  shear  wave  velocity  in  unconsolidated  mar¬ 
ine  sediments  are  sensitive  to  the  effects  of  compaction  and  the  decrease  of  void  ratio 
(Stoll  1989;  Richardson  et  al.  1991).  The  shear  strength  with  depth  also  increases 
from  1  to  12kPa  (Figures  5  and  6),  which  corresponds  well  to  shear  wave  velocity. 
Decrease  in  shear  strength  at  a  few  random  depths  may  be  artifacts  caused  by  the 
sampling  disturbance.  As  shown  in  Figures  3  and  4,  the  profiles  such  as  water  con¬ 
tent,  porosity,  density,  and  p-wave  velocity  do  not  reveal  significant  variation  with 
burial  depth,  compared  to  the  relative  variations  of  shear  properties.  It  is,  therefore, 
suggested  that  the  shear  properties  might  be  useful  to  predict  compaction  and  con¬ 
solidation  of  the  unconsolidated  sediments. 

Rigidity  modulus  shows  the  ranges  between  about  100  and  600dyn/cm2x  106 
(Figure  9).  The  relationship  between  rigidity  modulus  and  wet  bulk  density/porosity 
is  poor  (Figures  9a  and  9b),  whereas  rigidity  modulus  and  shear  strength/ 
shear  wave  velocity  shows  a  good  correlation  (Figure  9c,  9d)  because  shear  wave 
velocity  depend  mainly  on  grain  structure  and  moduli  of  the  sediment  frame 
(Hamilton  1971;  Orsi  and  Dunn  1991). 


Clay  Fabric  by  SEM 

Clay  fabric  (or  clay  microstructure)  is  defined  as  the  orientation  and  arrangement  or 
spatial  distribution  of  the  solid  particles  and  the  particle-to-particle  relationships 
(Bennett  et  al.  1991).  Fabric  changes  in  clay  are  mainly  related  to  consolidation,  min¬ 
eralogy,  and  grain  size,  as  well  as  diagenesis  (O’Brien  1970;  Bennett  et  al.  1981; 
Bryant  et  al.  1991;  Tribble  et  al.  1991).  The  microstructure  of  clay  fabric  strongly 
influences  and  largely  controls  the  physical  and  mechanical  properties  of  sediments, 
including  its  consolidation  behavior  (Bennett  et  al.  1991).  Clay  fabric  with  preferred 
orientation  provides  better  sediment  integrity  and  higher  shear  strength  because  of 
greater  surface  area  contacts,  and  higher  bonding  force,  compared  to  clay  sediment 
with  random  microstructure  that  has  lower  shear  strength.  Thus,  fabric  alteration 
appears  to  be  an  important  factor  influencing  both  shear  wave  velocity  and  shear 
strength  increase  of  the  core  sample.  This  relationship  between  shear  properties 
and  clay  fabric  is  observed  by  our  results  (Figures  5  to  8),  indicating  that  diagenesis 
may  play  a  significant  role  in  shear  properties  of  sediments  within  several  meters 
(Bennett  et  al.  1991).  The  microfabric  of  this  sediment  section  is  most  likely  in  the 
initial  stages  of  consolidation,  where  the  microstructure  has  sufficient  strength  to 
resist  the  stresses  of  the  small  overburden  load  impressed  upon  it  (Winterkorn 
1948;  Bolt  1956;  Lambe  1958;  Ingles  1968;  Bryant  et  al.  1991).  Accordingly,  the  cor¬ 
responding  pattern  in  the  variation  of  shear  properties  is  not  significantly  recognized 
from  the  geotechnical  properties  index  profiles  (Figures  3  and  4). 
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Rigidity  modulus,  (dynes/cm2*  106)  Rigidity  modulus  (dynes/cm2*  105) 


Figure  9.  Relationships  between  geotechnical  properties  and  acoustic  properties.  Linear 
regression  is  displayed.  Note  that  shear  properties  show  a  good  correlation. 


The  framboids  found  in  sediments  show  a  spheroidal  shape  (Figure  7b).  The 
spheroidicity  (spheroidal  morphology  of  framboids)  is  attributed  to  pseudomorph¬ 
ism  of  preexisting  spherical  bodies  (e.g.,  immiscible  organic  globules  or  infilling  of 
gaseous  vacuoles)  caused  by  bacteria  or  microflora  (Rickard  1970).  During  the  early 
diagenesis,  iron  diffusion  from  the  organic-rich  anoxic  sediment  into  the  bacterial 
cavity  resulted  in  iron  sulfide  precipitation  (Bennett  et  al.  1991).  The  authigenic 
pyrite  changes  the  clay  fabrics,  disturbing  the  particle  orientation  in  the  sediments. 
The  typical  framboid  found  commonly  in  organic-rich  sediments  displays  platy 
flakes  wrapped  around  the  sphere  and  is  associated  with  a  doming  of  adjacent  par¬ 
ticles,  as  shown  in  Figure  7b. 

The  samples  from  the  depths  of  50  cm  and  215  cm  with  different  values  of  shear 
properties  were  characterized  by  SEM,  as  shown  in  Figures  7  and  8.  For  example, 
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the  fabrics  of  the  samples  with  low  shear  strength  and  shear  wave  velocity  were  char¬ 
acterized  by  a  random  arrangement  of  clay,  while  the  clay  particles  of  sediments 
showing  high  shear  strength  and  shear  wave  velocity  were  well  oriented  caused  by 
increased  overburden  stress.  Therefore,  the  clay  fabric  study  can  be  a  major  tool 
in  understanding  the  mechanisms  of  the  consolidation  process  resulting 
from  porosity  reduction  and  its  relations  to  the  geotechnical  index  properties  of 
the  sediments. 

Conclusions 

Geotechnical  index  properties  (water  content,  porosity,  and  density)  with  the  excep¬ 
tion  of  upper  50  cm  and  below  200  cm  in  depth,  do  not  appear  to  decrease  with 
depth.  Compressional  wave  velocity  fluctuates  significantly  with  depth,  while  the 
shear  wave  velocity  (in  mean  value  of  horizontal  and  vertical  direction)  increased 
rapidly  from  8  to  20  m/s  in  upper  100  cm.  The  shear  strength  also  increased  from 
about  1  to  12kPa,  corresponding  to  the  shear  wave  velocity.  SEM  micrographs  of 
clay  at  215  cm  depth  show  the  preferred  particle  reorientation  compared  to  the  depth 
of  50  cm,  which  may  indicate  the  possible  processes  of  consolidation  affecting  the 
sediments.  Therefore,  the  consolidation  of  sediments  by  overburden  stress  causing 
the  clay  fabric  alteration  appears  to  play  an  important  role  in  changing  shear  proper¬ 
ties,  which  can  provide  important  information  to  predict  the  geotechnical  properties 
of  the  shallow  marine  sediments. 
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